Ground-state properties of a dilute Bose-Fermi mixture 
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We investigate the properties of a dilute boson-fermion mixture at zero temperature. The ground- 
state energy of the system is calculated to second order in the Bose-Fermi coupling constant. The 
Green's function formalism is applied to obtain results for the effect of Bose-Fermi interactions on 
the spectrum of phonon excitations, the momentum distribution, the condensate fraction and the 
superfluid density of the bosonic component. A quantitative discussion of these effects is presented 
in the case of mixtures with attractive interspecies interaction. 
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I. INTRODUCTION 

Recent experimental achievements in the trapping and 
cooling of mixtures of bosonic and fermionic atoms have 
attracted much attention on the physics of dilute Bose- 
Fermi mixtures in the quantum degenerate regime. Sta- 
ble Bose-Einstein condensates (BEG) immersed in a de- 
generate Fermi gas have been realized with ""Li in ^Li 
11 , 23Na in ^Li ||l and very recently with ^^Rb in '^°K 
In particular, the ®^Rb-^°K mixture appears to be 
highly interesting because of the large and negative in- 
terspecies scattering length which favours the stability 
of the mixture. Gurrent theoretical investigations of di- 
lute Bose-Fermi mixtures have mainly addressed the de- 
termination of the density profiles of the two compo- 
nents in trapped systems [Q, the problem of stability 
and phase separation and the BEC-induced inter- 
action between fermions H . The effect of boson-fermion 
interactions on the phonon excitation spectrum has been 
investigated in Refs. Q and Further, perturbation 
theory has been applied to the calculation of the ground- 
state energy in the regime of high fermion concentration 
by Albus et al. , whereas the opposite regime of low 
fermion concentration has been investigated long ago by 
Saam 0. 

In the present paper we consider a homogeneous di- 
lute mixture of a normal Fermi gas and a Bose-Einstein 
condensed gas at zero temperature. We calculate the 
ground-state energy of the system to second order in the 
Bose-Fermi coupling constant, interpolating between the 
high and low fermion concentration regimes, and obtain 
results for the dispersion of phonon excitations, including 
corrections to the sound velocity and Landau damping. 
In addition to these results, we focus on the effect of 
boson-fermion interactions on the momentum distribu- 
tion of the bosonic component, on the depletion of the 
condensate state and on the superfluid density. We find 
that the boson momentum distribution can be strongly 
affected by the fermions, resulting in a suppression of 
the occupation of low momentum states and in a long 
tail at large momenta. We also point out the possible 



occurrence of striking features caused by the Bose-Fermi 
coupling, such as the increase of the condensate fraction 
compared to the pure bosonic case and the superfluid 
density becoming smaller than the condensate fraction. 

The structure of the paper is as follows. In Sec. II we 
introduce the basic Hamiltonian of the system within the 
Bogoliubov approximation. In Sec. II A we calculate the 
ground-state energy by direct use of perturbation theory. 
In Sees. II B-D we introduce the Green's function for- 
malism which we apply to calculate the elementary exci- 
tation energies, the momentum distribution, the conden- 
sate fraction and the superfluid density of the bosonic 
component. In Sec. Ill we comment on the stability 
of the mixture against phase separation and we present 
numerical results obtained for specific configurations of 
mixtures with attractive Bose-Fermi coupling. Finally, 
in Sec. IV we draw our conclusions. 



II. THEORY 

We consider a stable homogeneous mixture of a Bose 
gas and a spin-polarized Fermi gas. The Hamiltonian of 
the system can be written as the sum 



H = Hp 



Hp 



Hi, 



(1) 



of the pure fermionic {Hp) and bosonic {Hb) Hamilto- 
nian, and of the interaction term {Hint) which accounts 
for the coupling between the two species. If one neglects 
p-wave interactions in the spin-polarized Fermi gas, Hp 
corresponds to the Hamiltonian of a free gas 



Hp 



k 



(2) 



where 6k, are the annihilation and creation operators 
of fermions and ef" = h^k^ /2mp is the energy of a free 
particle of mass mp. The pure bosonic component is 
described instead by the Bogoliubov Hamiltonian 



Hp 



Ep 



k 



(3) 
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written in terms of the quasi-particle annihilation and 
creation operators ak, aj^.. These operators are related 
to the bosonic particle operators Ok, a]^ through the well- 
known canonical transformation 



By using the above expression for , the total Hamilto- 
nian (|l|) takes the form 

H = Eb+ gBpNEUF + ^ Cfc &k^k + ^ t^fca^ak 



flk = MfcQ!k + Vka[ 

flk = Ukal. + WfcQ!_ 



(4) 



with coefficients uf, = 1 + v"^ — (e^ + 933^0 + LOk)/'^^k 
and UkVk = —gBBno/2ujk. The elementary excitation 
energies obey the usual Bogoliubov spectrum 



(5) 



where = fi^k^ /2mB is the free particle energy, ng is 
the condensate density and Qbb = ^Trh^a/mB is the cou- 
pling constant between bosons fixed by the mass tub and 
the s-wave scattering length a which we assume positive. 
The constant term 



AnnBa 
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2mBa 



(6) 
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The Hamiltonian (|lo| ) is treated in perturbation theory. 
We write H — Ho+Hint, where Hq — Hf+Hb is the un- 
perturbed Hamiltonian which is diagonal in the fcrmionic 
particle operators and in the bosonic elementary excita- 
tion operators, and 

Hint = QbeNbuf 

+ gBF\/»o y^(Mfc + vk){ai^ + a^k)P^k (11) 

k 

is the perturbation term. Since we make use of the pseu- 
dopotential approximation for both the boson-boson and 
the boson-fermion interatomic potential, the present ap- 
proach is valid if both gas parameters are small, i.e. if 
nsa^ < 1 and npb^ < 1. 



is the ground-state energy of a dilute Bose gas expressed 
in terms of the gas parameter n^a^, with ns = Nb/V 
the total density of bosons. Result (^) includes the zero- 
point motion of the elementary excitations. Finally, the 
Hamiltonian 



Htnt = gBF / dr nB(r)nF(r) 



(7) 



describes the interaction between fermions and bosons 
through the coupling constant gBF = ±27r?i^6/m/j, 
which can be either positive or negative depending on 
whether the interactions are repulsive or attractive. The 
coupling constant gBF is fixed by the reduced mass 
TOfl — mBiTiF / {mB + mp) and by the modulus of the 
boson-fermion s-wave scattering length b. By introduc- 
ing the bosonic and fcrmionic density fluctuation op- 
erators = 1/VV J dr e*'^ '"[nB(r) — and p^ — 
1/VV J dr e''' ''[ni?(r) — np], where uf = Nf/V is the 
mean density of fermions, the interaction Hamiltonian 
can be rewritten in the following form 



Hint = gBpNEUF + gBF ^ Pk P^h 



(8) 



The constant term in the above equation corresponds to 
the mean-field interaction energy, the second term de- 
scribes the coupling between density fluctuations in the 
two species. Within the Bogoliubov approximation the 
bosonic density fluctuation operator can be written as a 
linear combination of quasi-particle operators 



—y 



(ufe + Wfe)(ak + a^k) ■ (9) 



A. Ground-state energy 

The unperturbed ground-state energy can be readily 
calculated from the Hamiltonian Hq. One finds 



Eb + T^FeF , 
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(12) 



where Eb is given by (^) and €f — h'^kp/2mF — 
Tt? ['oTr'^nF)'^^^ /2mF is the Fermi energy, kp being the 
Fermi wave- vector. The mean- field term A£;(i) = 
gBF^BTiF provides the first order correction to E^. To 
second order in gBF one obtains the result 

k 



mHM\n)[^ 
Eo — En 



(13) 



where \n) are the excited states of the unperturbed 
system with energy En- The first term in the above 
equation arises from Ai?^^-* due to the renormalization 
of the boson-fermion scattering length gBF — *■ gBF + 
g^p\/V^^2mFi/'lT?k^ ■ The second term is the stan- 
dard result of second order perturbation theory. By di- 
rectly calculating the matrix elements of the interaction 
Hamiltonian Hint, one finds the following result for the 
ground-state energy 

E = Ea + NBgBpnF + NBg%Fi^Fj7 Y] 

k 

W 2 1 V-/ ^ ,2 <(l-<+k) 

^BgBFTH > + ^fe) 



k,q 



(14) 



|q+k| 
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where = d{kp — q) is the momentum distribution of 
the Fermi gas. We notice that the renormahzation of 
the boson-fermion scattering length is crucial in order 
to cancel the ultraviolet divergence in the integral over 
the wave-vector k. The dependence on the relevant pa- 
rameters becomes clearer by expressing the integrals over 
wave- vectors in units of kp. Result (|l^) takes the form 

E = Eq + NegBFriF + NBepiriFb^f/^Aiw, a) , (15) 

where the second-order energy shift is expressed in terms 
of the Fermi energy ep, the fermionic gas parameter npb^ 
and the dimensionless function 



A{w,a) = - 



3fc2(l + w) 



2/3 



1 



W 



dk 



+1 



(16) 



a 



dqq- 



1 - 61(1 - ^yY+W+2kqn) 

y/k^ + a + wk + 2qwQ 



In the above expression w = mB/mp is the mass ratio 
and a = 2/{kp£^B)'^ is a dimensionless parameter fixed by 
the product of the Fermi wave- vector kp and the Bose 
healing length £^b = 1/ y/^irriQa. In Fig. 1 we show the 
dependence of the function A{w, a) on the parameter a 
for a fixed value w = 1 of the mass ratio. Two regimes 
are worth studying at this point. The regime kp^B 3> 1 
(a ^ 1) corresponds to a system where the Fermi energy 
is much larger than the chemical potential of the bosons 
^ gBBiT-o (we assume w ~ 1). In this regime, expres- 
sion (^6[) can be simplified as it depends only on the mass 
ratio w: A{w,a — > 0) = Ai{w), where 



Ai{w) 



2 /eV^^ 1 + w 



3 Vtt 



dk / dn 

w Jo J-i 



(17) 



Jo 



(1 + iu)k + 2qwn 



In this regime result (^5|) reads 

E = Eo + NBgBFUp + NBep{npb^f'^Ai{w) , (18) 

and coincides with the finding of Ref. |jll| obtained us- 
ing the T-matrix approach. Notice that the effect of the 
Bose-Fermi interaction on the ground-state energy is in- 
dependent of the Bose-Bose interaction. In fact, in the 
limit kpS^B ^ 1, the relevant contribution to the inte- 
gral in Eq. (|lj) comes from the large momentum re- 
gion where the bosons behave as independent particles. 
The opposite regime kp^B <C 1 (a S> 1) corresponds 
to a Fermi energy ^ gBBno (w ~ 1). One finds: 
A(w, a — > oo) — 2^/aA2{w) , where 



A2iw) 



2 

3 Vtt 
1 



6^2/3^ 



w 



dk X 



il + w)k^ 



VTTk^{VTTW + wk) 



(19) 



In this regime the ground-state energy of the system can 
be written as 



E = Eo + NsgBFnp 



(20) 



+ ^T7t7^^-^-("-^^)^/^ (;^)'''(n,a3)V6A,(«;) . 



In contrast to result ( |18D , Bose-Bose interactions are im- 
portant if kp^^B ^ 1 and the second order correction de- 
pends explicitly on the Bose gas parameter. Result ( |20| ) 
has been first derived by Saam |l^ for a dilute mixture 
with a = b. The dependence of Ai and A2 on the mass 
ratio w is presented in Fig. 2. 
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FIG. 1. Dimensionless function A(w,a) [Eq. (p^J as a 
function of a for the value w; = 1 of the mass ratio. 




FIG. 2. Dimensionless functions Ai{w) [Eq. (p^j (solid 
line) and A2{w) [Eq. (|l|)] (dashed line). 
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B. Phonon spectrum 

In this subsection we calculate the effect of fermions on 
the dispersion of the bosonic phonon branch. To this aim 
it is convenient to use the Green's function formalism. 
The bosonic quasi-particle Green's function is defined as 



D(k,t) = --(r(akW4(0))), 



(21) 



where T(...) is the time ordered product. Perturbation 
theory provides us with a precise recipe for calculating 
the Green's function D{h, t) to a given order in Hint 



i 1 



(rK(i)4(o)^)> , 



(22) 



where the time evolution operator S is defined through 
the series expansion 



S = 



oo 

E 

n=0 



1 



+ 00 



dti.. 



dtr. 



X T 



Hint{il)---Hint {tn 



(23) 



In Fourier space the unperturbed quasi-particle Green's 
function is given by 



Co(k,^) 



1 



(24) 



where ry > and infinitesimally small. To second order 
in Hint a direct calculation gives the result 

Z?(k,c^) -I?o(k,c^) 

+ 9lFm{uk + VkfDl{]i,uj)Il^ (k,w) . (25) 

In the above equation 11^ (k, uo) is the density-density re- 
sponse function of an ideal Fermi gas defined as 



nak,c.) = i^n^(i 

q 



'q+k 



(26) 



1 



1q+k| 



flLO 



IT] 



By solving Eq. ( p5[ ) for D ^ (k, uj) one gets the following 
result for the excitation energy to order g^p 



Tiuj = LOk + g^p-ng— n^(k, huj 



(27) 



The real and imaginary part of the fcrmionic response 
function 11^ give rise respectively to a frequency shift 
and a damping of the quasi-particle. We are interested 
in the collective modes of the boson system which corre- 
spond to low momenta ujk ^ gBsno and k <^ kp. In the 
regime kpS^B S> 1 the fermionic density-density response 
function can be approximated by WM 



n^(k, fiUj ^ ujk) = 



mpkp 



l + i 



2V2wkF^B. 



(28) 



and one gets the following result for the phonon excita- 
tion energy 



huj = hkcB 



1 



(29) 



In the above expression the real part in the square 
bracket gives the correction to the Bogoliubov sound 
velocity cb — \/ gBBno/niB, and the imaginary part 
gives the Landau damping of the phonon excitation due 
to collisions with the fermions. We notice that the 
correction to the sound velocity (second term in the 
square bracket) can be rewritten as N{0)ggp/2gBB, with 
A^(0) — ■mpkp/{2Tr^h'^) being the density of states at the 
Fermi surface. The stability condition requiring that the 
real part of the excitation energy be positive, corresponds 
to the linear stability condition N{0)g^p/ gsB < 1 (see 
Sec. III). In the regime kF£,B ^ 1 the Bogoliubov 
phonon becomes soft when approaching the condition 
where the system phase separates. In the opposite regime 
^fC-B ^ 1, the imaginary part of 11^ vanishes at low 
momenta and the renormalized velocity of sound of the 
undamped phonon mode becomes 



huj — TikcB 



1 



^2 l2 



h np 



a? ub 



(30) 



The dispersion of phonons in dilute mixtures has also 
been discussed by Yip |^ under more general condi- 
tions using a numerical approach. He also predicts a 
shift towards lower frequencies with Landau damping for 
kpS,B 2> 1, and towards higher frequencies with no damp- 
ing for kpS^B < 1. 



C. Condensate fraction and momentum distribution 

The effect of the Bose-Fermi interaction on the momen- 
tum distribution of the bosons can be readily obtained 
by applying the perturbative approach to the bosonic 
single-particle Green's function 



Gn(k,t) = -^(T(ak(t)4(0))) 



(31) 



The unperturbed single-particle Green's function is given 
by 



G?i(k,c^) = uiDo{\^,uj) + vlD^{\^,-uj) 



huj — LUk + iv 'huj + LUk — 



(32) 



To second order in the Bose-Fermi coupling constant one 
gets the result 
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Gii(k,w) = G?i(k,w) +g|^no(Mfe +wfc)^n^(k,w) x 

X ulDl(k,Lu)+vlDl{k,-Lu)+ (33) 

+ 2ufcUfcL)o(k, w)L»o(k, -w) . 

The momentum distribution of the bosons can be ob- 
tained from the single-particle Green's function through 
the following relation 



ih 



which gives the result 



+00 



— / dc^Gii(k,c^)e'"\ (34) 



+ 



' ^|k+q| 

2ukVk 



0^ K + ' 



|k+q| 



1k+q| 



(35) 



By integrating over momenta one obtains the fraction 
of atoms which are scattered out of the condensate due 
to interaction. The quantum depletion of the condensate 
is given by 



N' 



Nb 30F 



(riija ) 



3U/2 



2 1 ef 1 
9bf77 2^—77 



(36) 



{uJk + 



k+q| 



cF\2 



|k+q| - '^q ) 



The first term in the above equation corresponds to the 
Bogoliubov depletion present also in a pure bosonic sys- 
tem which is due to interaction effects among the bosons. 
The second term accounts instead for boson-fermion scat- 
tering processes. The condensate fraction is obtained 
from the difference Nq/Nb = 1 — N' /Nb- By writing 
the integrals over wave- vectors in units of kp, result ( |3^ ) 
reads 

^ = -^{nBa')^^ + {npb'f/'Biw, a) . (37) 
Nb 3V7r 

As for result (p^), the dimensionless function B depends 
on the mass ratio w — niB/TUp and on the parameter 
a = 2/{kp^B)'^ and is given by 



B{w, a) 



12 /"^ f+ 



dqq- 



1 - 61(1 - ^q^ + k^ + 2kqn) 



\fWT~a Jo (Vfc^ + a + wk + 2qwVLY 

/ 1 a{wk + 2qwn)\ 

VVFT^" 2fc2(fc2 + a) ) ■ 



(38) 



The dependence of the function B on the parameter a is 
shown in Fig. 3 for the mass ratio w = 1. In the regime 
^fCs ^ 1 (a <C 1) the quantum depletion due to scat- 
tering processes with the fermions becomes independent 
of the Bose-Bose coupling constant and one finds 

^ = -^{uBa')'/' + {npby/^B,{w) , (39) 
Nb 3V7r 

where the function Bi of the mass ratio is given by 
12 



6l/37r2/3 



(w + l) 



dil / dqq 



(40) 



1 



{w + l)y/q^{n^ -l) + l + {w- i)qn 

In the opposite regime kpS^B ^ 1 (a ^ 1), the quan- 
tum depletion takes important contributions from the 
region of low momenta k < mcB/ft and interaction ef- 
fects among the bosons become crucial. In this regime 
one finds: B{w,a — > 00) = 2B2{w) / {3y/a) where the 
function B2 is given by 

12 



B2iw) = 



61/37^2/3 



(i + wY 



dk 



l + fc2 



1 



1 



2fcVl + fc2 



(41) 



k'^ + wfc)2 

and the quantum depletion takes consequently the form 
iV' _ 8 



3V^ 
TrVe b2 

62/3 ^2 



3U/2 



{nBa^)'/'^B2{w) 
ub 



(42) 



In Fig. 4 we plot the functions Bi and B2 of the mass 
ratio w. We notice that, while Bi is always positive, B2 
is positive for small values of w and becomes large and 
negative when ^ mp. Thus, if the fermions are 
light enough compared to the bosons, their effect in the 
regime kpS^B <C 1 is to enhance the occupation of the 
condensate state resulting in a decrease of the quantum 
depletion. Of course, when B2 is negative, result ( ^ ) is 
valid only if the second term is small compared to the 
first one. 
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FIG. 3. Dimensionless function B{w,a) [Eq. (pSj)] as a 
function of a for the value w = 1 of the mass ratio. 




FIG. 4. Dimensionless functions -Bi(to) [Eq. 
line) and B2(w) [Eq. @)] (dashed line). 
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Another quantity of interest is the momentum distribu- 
tion of the bosons given by eq. (|35|). For a pure system 
the momentum distribution is given by the well-known 
Bogoliubov result 



''k I pure 



+ 1 - x\/x^ + 2 
2xVx'^ + 2 



(43) 



which has been written in units of the inverse healing 
length X = kS^B- At small momenta x -^1, the momen- 
tum distribution is dominated by the quantum fluctua- 
tions of phonons which give rise to the infrared divergence 
Ipwe ^ l/\/8a;. At large momenta x S> 1, one finds 
instead the algebraic decay n^\pure l/4x'*. The impor- 
tant question arises concerning the effect of the fermionic 
component on . In units of the inverse healing length 
result (^ reads 



(67r2)i/3 

-1 rkp^B 

dfl / dyy^ X 

Vx^ + 2 J-i Jo 



^ 1 - eikpis - y/y' + x2 -I- 2xyn) ^ 
{^/x"^ + 2 ~\-wx + 2ywiiy 
1 wx + 2ywn\ 

V^^" x^x^ + 2) J • 



(44) 



In a mixture there are two characteristic wavenumbers 
in the problem: and kp. The structure of the dis- 
tribution is consequently richer and depends on how 
compares with kp- Nevertheless very interesting conclu- 
sions can be drawn in the two limits k <C (.g^^kp and 
k J> ^g^, kp. At low momenta {k <C ^B^jkp) the effect 



of the Bose-Fermi coupling is to enhance the occupation 
of the condensate by stimulating scattering of particles 
from k ^ states into the condensate. The contribu- 
tion to the momentum distribution is thus negative and 
reduces the coefficient of the 1 /x divergent term 



1 



1 



where 



VSx 
X f{2w/^) 



fit) 



61/3(1 + u;)2 5 

/i/3 -{npb^y/^ X 



1 + t t 



ln(l + t) 



(45) 



(46) 



Notice that the modulus of the Bose-Fermi contribution 
can be written as N{0)g^p/2gBB ■ fi'^w/^/a) [see com- 
ment after Eq. ([29|)]. Phase stability of the mixture re- 
quires that in any case N{0)ggp/gBB < 1- Since more- 
over < /(t) < 1 for all t > 0, we conclude that the 
Bose-Fermi term is always less, in modulus, than the 
Bose-Bose one and thus Eq. ( ^ ) is well defined. It is 
also worth noticing that in the regime kp^^B ^ 1 the 
function / ^ 1 and ^ (c/cb)/ V^x, where c is the 



renormalized sound velocity from Eq. (^9|). 

At large momenta (fc ^ ^B^jkp) the Bose-Fermi cou- 
pling provides instead further depletion of the condensate 
and one finds 



1 

4^ 



1-t- 



32 



3(6^2) 



2U/3 



[npb-') 



3N2/3"-a 



np 



(47) 



In the regime kp£,B 3> 1 the shape of can be very 
different compared to the pure case. In fact, the cut-off 
momentum increases from to kp and the momentum 
distribution acquires a long tail at large momenta (see 
Fig. 5). The effect of the Bose-Fermi interaction on the 
momentum distribution and condensate fraction of the 
bosons has been calculated using variational methods in 
the hypernetted-chain scheme for the strongly correlated 
liquid •^He - ''He mixture However, weak effects 

are found due to the low '^He concentration attainable in 
stable mixtures. 



D. Superfluid density 

An important aspect of Bose-Fermi mixtures concerns 
their superfluid behavior. In a pure bosonic system at 
zero temperature the density of the normal component 
vanishes, pn = 0, and the system is completely super- 
fluid. In a Bose-Fermi mixture the normal density does 
not vanish and can be written as p„ = rripUp in terms 
of the effective mass rup of the fermionic particles. In 
the absence of Bose-Fermi coupling, the effective mass 



coincides with the bare mass of the fermions, 



mp, 



and the normal density of the mixture coincides with 
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the mass density of the fermionic component. Interac- 
tion effects dress the fermions resulting in an effective 
mass larger than the bare mass. A well known exam- 
ple is provided by liquid '^He-'^He solutions, where due to 
the strong correlations the effective mass of '^He atoms 
is m*p/mF — 2.3 |Q. An equivalent way of interpreting 
the excess of normal density in Bose-Fermi mixtures is 
to write p„ = iriFnF + Pn, where is the contribution 
arising from the excitation of bosonic quasi-particles due 
to the Bose-Fermi coupling. In terms of the mass dif- 



ference Smp 



mp, the bosonic component of p„ 



is given by p^ — npSmp. The superfluid density of the 
system is defined as ps = pb ~ Pn ■ The bosonic compo- 
nent of the normal density can be calculated through the 
long- wavelength behavior of the static transverse current- 
current response function [p!5| 

Pn = — lim G'T(q, uj ~ 0) 



(48) 



In the above equation G'T(q, = —i/h{T{j^{t)j^^{0))) 
is the bosonic transverse current-current response func- 
tion, with q = (0, 0, q) directed along z and the x 
component of the current operator defined by = 

1 /y/V hkxol^ay^+q. In terms of single-particle Green's 
functions Gxifl,'-^) can be written as 



zvr V ^ — ' 

k 



2 / duj' X 



Gii{k,-uj')Gii{k + c[,uj - 
6*21 (k, a;')Gi2 (k + q,Lu - u 



(49) 



where G'i2(k,t) = ~i/h{T{ai,{t)a^k{0))) and G'2i(k,t) = 
— i/?i(T(aj^(t)aLj^(0))) are the anomalous single-particle 
Green's functions, which are peculiar of Bose condensed 
systems. The unperturbed anomalous single-particle 
Green's functions can be readily obtained and one finds 

GUk,u;)=GUKuj) 

= % ~~ ~ % — \ ~ ■ ^^^1 

TiLO — u}k+ irj nu! + uJk — "7 
To second order in ggp wc find the result 

Gi2(k,C^) =G2l(k,C^) 

= G?2(k,w) -|-g|^no(u/c + Vkfll^iKuj) x 
X {ul+vl)Doik,u;)DoiK~u;)+ (51) 

+ UkVkDl{k,uj) +UkVkDl{k,-LL>) . 



^From the perturbation expansions ( |33| ) and (^]J) one gets 
the following result for the long-wavelength limit of the 
static current-current response function 

lim GT(q,a; = 0) = ■7^j^y^^{fikx)'^g%pno{uk + VkY x 



27r V 

k 



4-00 



dujIi^{k,oj)Dl{k,uj) , (52) 



valid to order g\p. The fraction of Bose particles which 
contribute to the normal density of the system is thus 
given by the following expression 

PB r^'^v^ uok 



1^ <(i-<+k) 



-y 



(53) 



By expressing the integrals in units of the Fermi wave- 
vector kp one finds the result 



PB 



(54) 



in terms of the Fermi gas parameter and of the dimen- 
sionless function 



G(w7, a) 

fc2 



dqq 



1 - 61(1 - ^q^ -t- P 2kq^) 



\Jk'^ + a Jo (Vfc^ + a + wk + 2qwfl)^ 



(55) 



In Fig. 5 we show the dependence of the above function 
on the parameter a = 2/ (kp^B)'^ for the value w = 1 of 
the mass ratio. In the regime kp^B ^ 1 the function G 
becomes independent of the parameter a and result ( [5^ ) 
takes the form 



^ = (npb'f^'Ciiw) , 

PB 

where Gi is given by 

16 f^^ 

^l(^) = 61/3^2/3 dqq^ ^ 



(56) 



(57) 



{l + w)^q^{n^ -i) + i-qn 

[(1 + w)^q^n^ -l) + l + {w~ l)qfl\^ ■ 

In the opposite regime, kp^B <C 1, the function G 
takes the following asymptotic value C{w,a — s- oo) = 
2C2{w) / {3y/a) and one finds 



PB 6 ' a TiB 
where the function G2 is given by 

G^H^'-^^Tdk 



(58) 



61/3^2/3 VTTP(yrTfc2+wfc)3 



(59) 



Result (Hfl) is in agreement with the finding of Ref . px| ] 
where the fermion effective mass has been calculated for 
kp£,B <C 1 in the case of equal scattering lengths a = b. 
In Fig. 6 we show both Gi and G2 as a function of the 
mass ratio. 



7 




I — — — I 

0.01 0.1 1 10 

a 

FIG. 5. Dimensionless function C{w,a) [Eq. (^s])] as a 
function of a for the value ui = 1 of the mass ratio. 




0.1 1 10 



w 

FIG. 6. Dimensionless functions Ci{w) [Eq. dsst ] (solid 
line) and C2(m) [Eq. (H)] (dashed line). 



III. NUMERICAL RESULTS 

In deriving the above results we have been completely 
general. In this section instead, we look at their impli- 
cations on actual physical systems and we discuss more 
quantitatively the striking consequences of the presence 
of fermions on the Bose-Einstein condensate which we 
have pointed out. 

First of all we have to spend a few words on phase 
stability. The obvious requirement for all of the above 
to be valid is that bosons and fermions mix in the first 
place. According to the findings of ref. [|j one has to 
separately analyze two cases: attractive Bose-Fermi in- 
teraction, i.e. gsF < 0, and repulsive one, i.e. gsF > 0- 
The former case is simpler, as the only requirement is 
that arising from linear stability: N{0)ggp/gBB < 1, 
where iV(0) = mpkF /{'in^h'^) is the density of states at 
the Fermi surface. This condition fixes an upper limit on 
the fermionic density irrespective of the bosonic one 



477 /a\3 1 

If instead gsF is greater than zero, while the fermionic 
density cannot in any case be larger than nF,max, yet 
for every uf < riF.max there is a critical ns.dnF) at 
which the system phase separates anyway. The function 
nB,c{np) is non trivial, and one has to check specifically 
for a given pair of elements and densities whether the 
system is stable or not. 

For this reason it is certainly easier to observe impor- 
tant effects in a mixture where gs f is negative, as one can 
reach higher densities without running into problems. On 
the other hand, since our calculations are to order g^^p 
all effects computed here are identical for both signs of 

9BF- 

A mixture with negative gsF, and supposingly large 6, 
is already available experimentally. This is the ''"K-^^Rb 
mixture recently reported in ref. [^. According to the 
authors h — 300(100)ao in units of the Bohr radius. 

We consider a mixture of ^^Rb and ''°K with densities 
ub = nF — 3 ■ 10^^ cm^'^. For the scattering lengths 
we use a =110ao and b =300ao, and the mass ratio is 
w =2.175. The value of the parameter kp^B is 3.9. 

In Fig. 7 we show results for the boson momentum 
distribution of the mixture, obtained from Eq. ( ^ ) with 
the parameters given above, compared with the one of a 
pure system [eq. (^)]. The total value of the quantum 
depletion is ~ 3% for this configuration and the Bose- 
Bose and Bose-Fermi contributions are ~ 1% and ~ 2% 
respectively (see Fig. 8). We see a large effect due to 
boson- fermion interactions. The momentum distribution 
at low momenta is depressed while at large momenta is 
considerably enhanced and a tail appears which extends 
to very large momenta. 

In Fig. 8 we report results on the dependence of the 
condensate depletion and normal fluid fraction as a func- 
tion of the fermion density np- The boson density has 
been fixed to 3 • 10"'^'' cm"'^. With the above values of the 
parameters the system becomes unstable at a fermion 
density of 5.2 • 10^^ cm~^. We see that, by increasing the 
fermion density, both the condensate depletion and the 
normal fluid fraction increase. The prediction for fermion 
densities close to the instability threshold gives N' /Nb ~ 
4% and p^/pB 2%. 

In Fig. 9 we show the condensate depletion and the 
normal fluid fraction of a hypothetical mixture with mass 
ratio rriB/'mp = 0.1. We also assume gBF < and 
we have fixed the scattering lengths to a = IOOoq and 
b — 50ao, and the boson density to ub — 3 ■ 10^^ cm^'^. 
Such system becomes unstable for np > 1.0 • 10^^ cm~'^. 
We see that an increase oinp is followed by an increase of 
both the condensate depletion and the normal fluid frac- 
tion. However, at the fermion density np ~ 2-10"'^^ cm~'^, 
the normal fluid fraction becomes larger than the conden- 
sate depletion. This effect has been first discussed in con- 
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nection with disordered dilute Bose gases by Huang and 
Meng [ [T6t and very recently has been explicitly proved 
in a Bose gas with quenched impurities using quantum 
Monte-Carlo techniques It is striking that the same 
effect is found for Bose- Fermi mixtures in regimes of high 
fermion concentration when mp ^ . 

Finally, in Fig. 10 we present results of the condensate 
depletion and normal fluid fraction of a different hypo- 
thetical mixture with mass ratio ms/mj? = 20. We as- 
sume again an attractive Bose- Fermi coupling {gsF < 0) 
and we use the following parameters: a = lOOao, b = 
150ao and = 3 • 10^* cm^'^. The mixture becomes 
unstable for np > 2.3 • 10^^ cm~'^. In this case we see 
a different striking behavior: for low fermion concentra- 
tion the effect of the Bose- Fermi coupling is to stimulate 
the occupation of the condensate, resulting in a quantum 
depletion smaller than in the pure case. 

Due to the smallness of the effects, the interest of the 
results reported in Figs. || and |l^ is to show that they 
can actually occur in realistic mixtures, though it might 
be difficult to observe them with present techniques. 




FIG. 7. Boson momentum distribution of a *^Rb -^"K mix- 
ture (solid line). The momentum distribution of a pure system 
(short dashed) and the Bose- Fermi contribution (long dashed) 
are also shown. In the figure as in the text x stands for fc^s. 



IV. CONCLUSIONS 

We have studied the properties of a dilute Bose-Fermi 
mixture at zero temperature using a perturbation ap- 
proach. We have investigated both equilibrium prop- 
erties, such as the ground-state energy, the boson mo- 
mentum distribution and the normal fluid fraction, and 
dynamic properties such as the dispersion law and damp- 
ing of phonon excitations. The system is very rich. By 
varying the mass ratio, the densities and the scattering 
lengths of the two components one can obtain very differ- 
ent regimes where striking effects due to boson-fermion 
interactions can occurr. These include: the strong sup- 
pression of the boson momentum distribution at low mo- 
menta, localization effects exhibited by the superfluid 
density which becomes smaller than the condensate frac- 
tion and stimulated scattering of bosons into the conden- 
sate. 




le+13 



FIG. 8. Condensate depletion and normal fluid fraction of 
a '^'^Rb-'^^K mixture as a function of the fermion density. The 
total condensate depletion is shown with the solid line. The 
Bogoliubov and the fermion depletion are shown respectively 
with the horizontal short-dashed line and the long-dashed 
line. The dotted line corresponds instead to the normal fluid 
fraction. 
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FIG. 9. Condensate depletion (solid line) and normal fluid 
fraction (dotted line) of a hypothetical mixture with mass 
ratio ms/mi? =0.1 as a function of the fermion density. 
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FIG. 10. Condensate depletion and normal fluid fraction 
of a hypothetical mixture with mass ratio mB/mp =20 as a 
function of the fermion density. The line codes used are the 
same as for Fig. 8. 
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